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1. Introduction 

Asymmetric systems offer unique information about the underlying dynamics, not available 
in symmetric proton-proton or nucleus-nucleus collisions. The asymmetry manifests itself in an 
asymmetric distribution of charged particles with respect to zero rapidity (or pseudorapidity) as 
measured by BRAHMS ^ and PHOBOS The asymmetry can be quantified by introduc- 
ing the ratio of pseudorapidity densities at a given negative pseudorapidity relative to that at the 
positive pseudorapidity of the same magnitude. This backward/forward ratio (forward being the 
original direction of motion of the light partner) is referred to as pseudorapidity asymmetry. The 
STAR Collaboration published pseudorapidity asymmetries in 200 AGeV dAu collisions for sev- 
eral identified hadron species and total charged hadrons in the pseudorapidity intervals |t]| < 0.5 
and 0.5 < < 1.0 [^. Asymmetries with the backward/forward ratio above unity for transverse 
momenta up to w 5 GeV/c are observed for charged pion, proton-i-anti-proton, and total charged 
hadron production in both rapidity regions. We anticipate that proton-lead (or deuteron-lead) data 
will soon be collected at higher energies, at the Large Hadron Collider (LHC). 

In the present study we investigate the roles of nuclear shadowing and multiple scattering 
in the generation of nuclear modifications and pseudorapidity asymmetries in a wide transverse- 
momentum range. We use the HUING shadowing parameterization [Q] and the Eskola-Paukkunen- 
Salgado (EPS08) nuclear parton distribution functions (nPDFs) While the former has been 
applied widely, the latter was not available at the time of similar earlier studies. We present pseu- 
dorapidity asymmetries for pBe at 30.7 GeV (Fermilab) and dAu at 200 AGeV (RHIC), where data 
are available from the E706 experiment [^, PHENIX [Q], and STAR [^]. We then make predictions 
for dPb at 8.8 ATeV (LHC). 



2. Model framework 

The invariant cross section for the production of final hadron h from the collision of nucleus A 
and nucleus B {A + B^h + X), can be written as 

Eh—^ = £ d^bd^rtA{b)tB{\b-r\)dXadXbdkTadkThdZcfa/A{Xa,kTa,Q^) 
^ P abed-' 

fh/B[Xb,ICTb,Q ) :^—^s5[S + t + U) , (2.1) 

' dt nzi 

where all quantities have their usual meaning [^]. 

The coUinear parton distribution functions (PDFs) are generalized to include a transverse mo- 
mentum degree of freedom, kj, as required by the uncertainty principle. This can be formally 
implemented in terms of unintegrated PDFs [§, |T^. To avoid some of the complications asso- 
ciated with using unintegrated PDFs, we resort to a simple factorized approximation, where the 
^r-broadened parton distribution in the nucleon is written as 

faMx,kr,Q') — 8{kT)-faMx,Q\ 8{kT) = '^P^-fT/ikhpp) ^ (2.2) 
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with f„ij^ (x, ) denoting the standard collinear PDF in the nucleon, and {k^ )pp is the two-dimensional 
width of the transverse-momentum distribution in the proton, related to the magnitude of the aver- 
age transverse momentum of a parton by {kj)pp = 4{kT)^p/7i. 

Most shadowing parameterizations include at least some of the effects of multiple scattering 
in the nuclear medium, while the HIJING parameterization (as we discuss further in Sec. ^ needs 
to be augmented with modeling nuclear multiscattering. For this purpose, in pA collisions we use 
a broadening of the width of the transverse momentum distribution according to 



{l4)pA = {kj)pp+C hpA{b) 



(2.3) 



where {kj)pp is the width already present in proton-proton collisions, hpA{b) is the number of 
effective nucleon-nucleon {NN) collisions as a function of nucleon impact parameter b, and C is 
the average increase in width per A^A^ collision. 

The collinear nPDFs fu/A{^:Q^) expressible as convolutions of nucleonic parton distri- 
bution functions (PDFs) fu/Ni^^Q^) a shadowing function ^'^/^(x, 2^) which encodes the nu- 
clear modifications of parton distributions. We use the MRST2001 next-to-leading order (NLO) 



PDFs [ ]1 1[ ] for the nucleon parton distributions and for the shadowing function we employ both 
the EPS08 shadowing routine [Bl] and HIJING [|]. (Other nPDFs, like FGS HKN [~ 



and 



the earlier EKS |14], are used elsewhere to calculate pseudorapidity asymmetries [15].) For the 



final hadron fragmentation we utilize the fragmentation functions in the AKK set 016|]. The factor- 



ization and fragmentation scales are from the best-fit results obtained in Ref. Jl7|]. We obtain the 
density distribution of the deuteron from the Hulthen wave function [111] (as in Ref. [19]), while a 



Woods-Saxon density distribution is used for gold and lead with parameters from Ref. [20]. 



3. Nuclear modifications and pseudorapidity asymmetry 

The nuclear modification factor, R\g{pT,'f]), and pseudorapidity asymmetry, Y^^^,^^{pt), can 
be defined for any produced hadron species h as 

1 Ehd^a^s/d^p\^ 



RUPT,ri] 



{Nbin) Ei4'a';p/d^p\r^ 



£j3 

dp 



Eh 



d^CT^ 



d^p 



(3.1) 



where {Nhin) is the average number of binary collisions. 

Let us consider the (double) ratio of the forward and backward nuclear modification factors in 
dAu collisions for species h: 



' RdAuiPT,ri) Ehd'a<^p/d^p\-^ / Ehd'a';p/d^p\r, 



(3.2) 



As discussed in Ref. [pT]], the pp rapidity distribution is symmetric around j = and thus cancels. 
We therefore obtain: 

RliAuiPT-r]) 



^AsymiPT) - R\{Pt) 



RdAuiPT,n) 



(3.3) 
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4. Nuclear shadowing and multiple scattering 

We expect that particle production in pA and dA collisions will be different in the forward and 
backward directions. This is because the respective partons have different momentum fractions 
(shadowing differences) and because the forward-going parton has to traverse a large amount of 
matter. In this Section we therefore study the interplay of shadowing and multiple scattering in 
more detail using the HIJING parameterization. 




Figure 1: (Color Online) Pseudorapidity distribution for minimum bias neutral pion production from dAu 
collisions at 200 AGeV. The solid line represents HIJING shadowing with intrinsic fcj in the proton, the 
dashed line is HIJING with intrinsic kj and multiple scattering, while the dot-dashed curve is obtained by 
turning off any transverse momentum. 

In Fig. [l] we display the pseudorapidity distribution for neutral pion production from dAu col- 



lisions at a fixed pT (= 2.25 GeV/c). The parameters of eq. ( |2.3| ) are unchanged from our previous 
studies at midrapidity [|^, [I^, and we have chosen a transverse momentum value comparable to 
(^t)pp, where the various effects are clearly displayed. 

The results of our study can be summarised as follows: In the absence of shadowing, the dis- 
tribution is symmetric around midrapidity. With shadowing only, an asymmetry develops: the yield 
at a fixed negative pseudorapidity, say T] = —2 {Au side or "backward") is higher than at the corre- 
sponding positive pseudorapidity (y = 2,d side, or "forward"). Thus the pseudorapidity asymmetry, 
YAsym, is greater than unity in this case. The inclusion of the intrinsic transverse momentum in the 
proton, while increasing the yields on both sides, does not affect the forward/backward symmetry. 
When multiple scattering is included, the yield shows a stronger increase in the forward direction, 
leading to YAsym less than unity. Thus we can say that shadowing suppresses the yield more on 
the d side (forward) relative to a symmetric collision, while the multiple scattering contribution is 
understandably large in the forward direction. 

We have carried out similar studies at higher pj values. When kj << pr, then kj effects 
become naturally smaller. Shadowing effects also become smaller as the antishadowing region 
of the HIJING parameterization is approached. Thus, at > 15 GeV/c the influence of multiple 
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scattering and intrinsic kj become negligible. These phenomena are most important at intermediate 
Pt values, 2 GeV/c ^ ^ 8 GeV/c at RHIC. It is interesting to note that a similar transverse- 
momentum region is sensitive to nuclear effects at lower (e.g. CERN SPS) energies, due to the 
~ log(y(s) scaling of the Cronin peak [23, |1, 25]. 



5. Pseudorapidity asymmetry 

5.1 Asymmetry in pBe collisions at 30.7 GeV 

The pseudorapidity asymmetry for the rapidity interval 0.2 < |t]| < 0.7 is shown in the upper 
panel of Fig. | compared with the E706 data Both EPS08 with kj and HIJING without multiple 
scattering give very small asymmetries and the data are also consistent with YAsym = 1 for low 
Pt- The HIJING parameterization with multiple scattering yields somewhat lower values at all 
transverse momenta. This is due to multiple scattering moderately increasing the yield on the /j-side 
relative to that of the Be-side. In view of the rather large error bars, all three sets are in reasonable 
agreement with the data. The lower panel is our prediction for the interval 1.0 < |t]| < 1.5. The 
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Figure 2: (Color Online) Pseudorapidity asymmeti-y, YAsym for +Be ^ tt" +X at 0.2 < \r] \ < 0.7 (top) 
and 1 .0 < \rj \ < 1 .5 (bottom). The solid line represents the EPS08 nPDFs, while the dashed line is obtained 
from HIJING with the inclusion of multiscattering. The dot-dashed line corresponds to HIJING without 
multiscattering, and filled triangles denote the E706 data 1^]. 

calculated effects are larger, but with similar structure to those at lower rj . 

5.2 Asymmetry in dAu collisions at 200 A GeV 

Figure ^ shows the pseudorapidity asymmetry for tt" production from dAu collisions at RHIC, 
for different pseudorapidity intervals. The two uppermost panels are our results for the asymmetry 
at |t]| < 0.5 and 0.5 < |t]| < 1.0 compared with the STAR data For pr > 4.0 GeV/c, the 
agreement with data is quite good for all three sets. At lower pr, multiple scattering increases 
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the calculated yield mostly in the forward direction as discussed in Sec. leading to calculated 
asymmetries below unity. At very high pT we observe a divergence in the model predictions. The 
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Figure 3: (Color Online) Pseudorapidity asymmetry, Y^sym f^i" d +Au ^ Tt^ +X different pseudorapid- 
ity intervals. The solid line represents the EPS08 nPDFs, while the dashed line is obtained using HIJING 
shadowing with the inclusion of multiscattering. The dot-dashed line corresponds to HIJING without multi- 
scattering, and filled triangles denote the STAR data 

lower four panels are our predictions for the asymmetry as pseudorapidity increases. The first three 
of these correspond to the BRAHMS pseudorapidity intervals [^6|]. The general trend is that the 
asymmetry becomes larger as T] increases. This mainly arises from the sti^ong shadowing in the 
larger nucleus at lower x values. Also, since increasing T] leads to decreasing accessible pt due to 
phase space constraints, the effects of multiple scattering become more pronounced. 



5.3 Asymmetry in dPb collisions at 8.8 ATeV 

Let us now turn to our predictions for the pseudorapidity asymmetry at the LHC energy of 
8.8 ATeV. The calculated results are displayed in Fig. ^ where the upper panel is for the interval 
ItjI < 0.9 and the lower panel is for 2.4 < ItjI < 4.0. These intervals correspond to acceptance in 



the central detector and in the muon arm, respectively, of the ALICE experiment [|27p. All three 
sets predict minimal asymmetry of the order of a few percent for the interval |tj| < 0.9. As we 
move to higher T] , the predicted asymmetry becomes more significant. As can be seen in the lower 
panel of Fig. ^, both EPS08 and HUING predict substantial asymmetry. According to the model, 
an LHC measurement of pseudorapidity asymmetry at around pj ~ 40 GeV/c could discriminate 
between the EPS08 and HHING shadowing prescriptions. 

At the present level, neither model variant gives agreement with all aspects of the data: in 
an earlier calculation we have found that shadowing parameterizations which do not need to be 



augmented by a multiple scattering prescription [14, 13, 12] have difficulty describing central 
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Figure 4: (Color Online) Predicted pseudorapidity asymmetry, YAsym for d + Pb +X at ^/s = 8.8 

ATeV for I?] I < 0.9 and 2A<\r]\< 4.0. The solid line represents the EPS08 nPDFs, while the dashed Hne 
is obtained from HIJING with the inclusion of multiscattering. The dot-dashed line corresponds to HIJING 
without multiscattering. 



to-peripheral ratios at forward rapidity 1 15]. We have checked that this also holds for the EPS08 
nPDFs. On the other hand, the HIJING parameterization with multiscattering yields pseudorapidity 
asymmetries below unity at low transverse momenta. 



6. Conclusion 

The mechanisms for particle production in asymmetric collisions leads to observable asym- 
metries in the pseudorapidity distribution of the produced particles at some collision energies and 
transverse momenta. We have considered the effects of nuclear shadowing and multiple scattering 
on pseudorapidity asymmetry for three asymmetric systems: pBe, dAu, and dPb in a wide cm. 
energy range from FERMILAB up to LHC energies. 

Overall, the calculated asymmetries are in reasonable agreement with available experimental 
data. Intrinsic transverse momentum in the nucleon is seen to be important at low pj. Multiple 
scattering increases the yield in the "forward" or positive pseudorapidity region, thus leading to 
a tendency for asymmetries less than unity at low pj in a scheme explicitly relying on multiple 
scattering, at variance with the data. An LHC measurement of pseudorapidity asymmetry in dPb 
collisions may be able to distinguish between the EPS08 and HIJING nPDFs at high pseudorapidity 
and high pj. 

A major constraint in assessing pseudorapidity asymmetries is the limited availability of data 
for direct comparison with theoretical calculations. More data in asymmetric light-on-heavy colli- 
sions separated with respect to positive and negative pseudorapidities are needed to judge calculated 
pseudorapidity asymmetries. At RHIC, the data from the high-statistic dAu Run 8 can be used to 
provide such a large data sample. 



7 



Pseudorapidity asymmetry in p(d)Pb collisions at the LUC 



George Fai 



7. Acknowledgments 

This work was supported in part by Hungarian OTKA PD73596, T047050, NK62044, and 
IN71374, by the U.S. Department of Energy under grant U.S. DOE DE-FG02-86ER40251, and 
jointly by the U.S. and Hungary under MTA-NSF-OTKA OISE-0435701. 

References 

[1] I. Arsene etal. [BRAHMS Collaboration], Phys. Rev. Lett. 94, 032301 (2005). 

[2] B. B. Back et al. [PHOBOS Collaboration], Phys. Rev. C 72, 031901 (2005). 

[3] B. 1. Abelev et al. [STAR Collaboration], Phys. Rev. C 76, 054903 (2007); J. Adams etal. [STAR 
Collaboration], Phys. Rev. C 70, 064907 (2004). 

[4] S. J. Li and X. N. Wang, Phys. Lett. B 527, 85 (2002). 

[5] K. J. Eskola, H. Paukkunen and C. A. Salgado, arXiv:0802.0139 [hep-ph]. 

[6] L. Apanasevich et al. [Fermilab £706 Collaboration], Phys. Rev D 68, 052001 (2003). 

[7] S.S. Adler et al. [PHENIX Collaboration], Phys. Rev. Lett. 91, 072303 (2003); Phys. Rev. C74, 
024904 (2006). 

[8] Y. Zhang, G. I. Fai, G. Papp, G. G. Bamafoldi and P Levai, Phys. Rev. C 65, 034903 (2002). 
[9] J. C. Colhns, T. C. Rogers and A. M. Stasto, Phys. Rev. D 77, 085009 (2008). 
[10] M. Czech and A. Szczurek, J. Phys. G 32, 1253 (2006). 

[11] A. D. Martin, R. G. Roberts, W. J. Stirhng and R. S. Thome, Eur. Phys. J. C 23, 73 (2002). 
[12] L. Frankfurt, V. Guzey and M. Strikman, Phys. Rev. D 71, 054001 (2005). 

[13] M. Hirai, S. Kumano and T. H. Nagai, Phys. Rev. C 70, 044905 (2004); Nucl. Phys. Proc. Suppl. 139, 

21 (2005). 

[14] K. J. Eskola, V. J. Kolhinen and C. A. Salgado, Eur. Phys. J. C 9, 61 (1999). 

[15] A. Adeluyi and G. Fai, Phys. Rev C 76, 054904 (2007). 

[16] S. Albino, B. A. Kniehl and G. Kramer, Nucl. Phys. B 725, 181 (2005). 

[17] P Levai, G. G. Barnafoldi, G. Fai and G. Papp, Nucl. Phys. A 783, 101 (2007). 

[18] L. Hulthen and M. Sugawara, "Handbuch der Physik", vol. 39 (1957). 

[19] D. Kharzeev, E. Levin and M. Nardi, Nucl. Phys. A 730, 448 (2004) [Erratum-ibid. A 743, 329 
(2004)]. 

[20] C. W. De Jager, H. De Vries and C. De Vries, Atom. Data Nucl. Data Tabl. 14, 479 (1974). 

[21] G. G. Bamafoldi, A. Adeluyi, G. Fai, P Levai and G. Papp, arXiv:0807.3384 [hep-ph]. 

[22] P Levai, G. Papp, G. G. Bamafoldi and G. I. Fai, Eur. Phys. J. ST 155, 89 (2008). 

[23] L. Apanasevich et al, Phys. Rev. D 59, 074007 (1999); L. Apanasevich etal. [Fermilab E706 
Collaboration], Phys. Rev. Lett. 81, 2642 (1998). 

[24] M. Zielinski, arXiv;hep-ph/981 1278. 

[25] G. G. Bamafoldi, P Levai, G. Fai, G. Papp and B. A. Cole, Int. J. Mod. Phys. E 16, 1923 (2007). 

[26] 1. Arsene et al. [BRAHMS Collaboration], Phys. Rev Lett. 93, 242303 (2004). 

[27] B. Alessandro et al. [ALICE Collaboration], J. Phys. G 32, 1295 (2006). 



8 



